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Abstract 

We study a parametric class of isotropic but not necessarily stationary Poisson hyperplane 
tessellations in n-dimensional Euclidean space. Our focus is on the volume of the zero 
cell, i.e. the cell containing the origin. As a main result, we obtain an explicit formula for 
the variance of the volume of the zero cell in arbitrary dimensions. From this formula we 
deduce the asymptotic behaviour of the volume of the zero cell as the dimension goes to 
infinity. 
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1 Introduction 

The majority of contributions to random tessellations is devoted to investigations in low and 
fixed dimensions. In particular, there exist only a few results on random tessellations in high 
dimensions, that is, with focus on asymptotic aspects as the dimension goes to infinity. Recently, 
the typical cell of a stationary Poisson- Voronoi tessellation in high dimensions has been studied 
in p], [26] and |37J. Alishahi and Sharifitabar [1] investigate the asymptotic behaviour of the 
volume and the shape of the typical cell of a stationary Poisson- Voronoi tessellation as the 
dimension n of the space goes to infinity. In particular, they showed that the variance of 
the volume of the typical cell converges to zero exponentially fast as n — > oo whereas it is 
well known that the expected volume is independent of the dimension. In the course of their 
investigation, they made use of an explicit formula for the variance of the volume of the typical 
cell in arbitrary dimensions. The asymptotic behaviour of the volume of the typical cell was 
studied earlier in the more general context of the nearest neighbour analysis by Newman et 
al. in [2H] and [2H]. They already showed convergence in distribution, though not via the 
convergence of the variance. 

In this work, we consider a parametric class of Poisson hyperplane tessellations and focus 
on the volume of the cell containing the origin (the zero cell). It is then natural to explore 
whether an asymptotic behaviour similar to that of the typical cell of a stationary Poisson- 
Voronoi tessellation is exhibited by the zero cell in the present more general class of random 
tessellations. An interesting family of not necessarily stationary or isotropic Poisson hyperplane 
tessellations is introduced in connection with the investigation of Kendall's conjecture in [10J. 
In the isotropic case, these Poisson hyperplane tessellations are completely determined by two 
parameters, the intensity 7 G (0, 00) and the distance exponent r e (0, 00). For a special 
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choice of the distance exponent and the intensity, the zero cell is equal in distribution to the 
typical cell of a stationary Poisson-Voronoi tessellation (cf. pU], [U Sections 5.2 and 5.3.2]). 
Therefore this family of isotropic Poisson hyperplane tessellations provides a general framework 
for investigating tessellations in high dimensions. Trying to extend the approach of Alishahi 
and Sharifitabar pQ to the wider context of the zero cell of this family of Poisson hyperplane 
tessellations, we came across the need for a formula for the variance of the volume of the zero 
cell. Finding a manageable expression turned out to be a rather complex issue. In fact, the 
formulas presented here mark the starting point of a more detailed study of the asymptotic 
behaviour of characteristics of the zero cell as the dimension goes to infinity. Results concerning 
lower dimensional sections of the zero cell and other shape characteristics as well as a connection 
to the hyperplane conjecture will be considered separately. 

In the following, we give a more detailed overview of our results. A precise description of the 
particular (parametric) model of a Poisson hyperplane tessellation used here is given in Section 
2. For this model we then derive, in Section 3, an explicit expression for the expectation 
and estimates of the moments of the volume of the zero cell in Proposition [T] An explicit 
expression for the second moment and the variance of the volume of the zero cell is provided 
in Theorem [TJ These results follow from a sequence of lemmas which make use of integral 
geometric transformations and the symmetries of the geometric situation. In Theorem [2] we 
deduce estimates for the variance of the volume of the zero cell which involve auxiliary quantities 
D(n, r) and E(n, r). These quantities have to be evaluated and estimated for a given distance 
exponent r and an intensity 7. In Corollary [T] we consider the zero cell of a Poisson hyperplane 
tessellation with constant distance exponent r. The choice r = 1 corresponds to a stationary 
Poisson hyperplane process. For constant intensity we then prove that all moments as well 
as the variance of the volume of the zero cell converge to infinity as the dimension n goes to 
infinity. 

In order to fix the expected volume of the zero cell, independent of the dimension, the 
intensity of the underlying Poisson hyperplane process can be chosen appropriately as a function 
of the dimension n. However, it follows from our estimates that as long as the distance exponent 
r is fixed, the variance of the volume of the zero cell is still divergent as n goes to infinity. The 
investigation in Section 4 thus suggests that in order to ensure that the variance converges 
to zero, the distance exponent r has to be adjusted to the dimension n. In Corollary [2] we 
summarize the case where the distance exponent r is proportional to the dimension n, i.e. 
r = an with a fixed factor a G (0, 00). For constant intensity we show that the expectation and 
the moments of the volume of the zero cell now all converge to zero as the dimension n goes to 
infinity. 

In Theorem [3] we finally consider the situation where the distance exponent r is proportional 
to the dimension, i.e. r = an with some fixed a > 0, and the intensity 7(0,71) is chosen as a 
function of the dimension n and the factor a in such a way that the expected volume of the 
zero cell is equal to a positive constant. In this case we prove that the variance of the volume 
of the zero cell converges to zero at an exponential speed as n — > 00. In particular, the 
volume of the zero cell converges in distribution. In the special case r = n (i.e. a = 1), we 
recover results for the typical cell of a Poisson-Voronoi tessellation. Of course, these findings 
are consistent with the results obtained in [Tj, though the estimates for the variance found by 
Alishahi and Sharifitabar are sharp. The present more general approach applies to a larger 
class of tessellations and admits various other extensions and variations that will be discussed 
in detail in subsequent work. 
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2 Preliminaries 



In the following, we mainly use the notation and terminology of the monograph [M]. We work 
in an n-dimensional real Euclidean vector space IR n , n > 2, with scalar product (•, •) and norm 
|| ■ || . The unit ball {x G lR n : \\x\\ < 1} centred at the origin o is denoted by B n , its boundary is 
the unit sphere S^ 1 . For k G {0, . . . , n}, the Grassmannian of /c-dimensional linear subspaces 
of IR n is denoted by G(n,k), and the afflne Grassmannian of /c-dimensional affine subspaces 
(Axflats) by A(n,k); both are equipped with their standard topologies. For u G S 11 ^ 1 and 
t G [0, oo), we write 

H(u,t) := {x G R n : (x,u) = t}, H~(u,t) := {x G E" : (x,u) < t}. 

Lebesgue measure on M n is denoted by A. For E G G(n, k), Lebesgue measure on E is denoted 
by \ E . Besides we define S^ 1 := E n S"^ 1 and H E (u,t) := E n H(u,t), for u G S 1 ^' 1 and 
k G {l,...,n}. The s-dimensional Hausdorff measure is denoted by "H s , where s > 0. For 
s = n we sometimes refer to it as n-dimensional volume V n . A frequently occurring constant is 
the volume of the unit ball, 

The surface area of the unit sphere S n ~ l is given by 

2vr? 

Un ■=n n -\S n - 1 ) =nn. 



For m G N and Xi, . . . ,x m G M. n , we denote by [xi, . . . , x m ] the convex hull and by span{xi, . . . , x„ 
the linear hull of 

The family of nonempty, compact, convex subsets of W 1 is denoted by /C n . For a topological 
space (T, T), a measure is always defined on the a-algebra £>(T) of Borel sets of T, i.e. the 
smallest a-algebra containing the open sets T. We write fi k := /i <E> ■ ■ • <8> /x, with factors /x, 
for the /c-fold product of a measure /x. By SO n we denote the group of proper rotations on M. n , 
and u n is the unique Haar probability measure on SO n . 

The following setting has previously been considered in a more general, not necessarily 
isotropic framework, in the context of Kendall's problem [10] (see also [6], jl]). Let (Q,A,P) 
denote the underlying probability space. Further, let X be a Poisson hyperplane process in 
M. n , i.e. a Poisson point process in the space A(n,n — 1). Subsequently, we identify a simple 
counting measure with its support, so that for a Borel set A C A(n,n — 1) both notations 
X(A) and card(XflA) denote the number of elements of X in A. We assume that the intensity 
measure O(-) = EA(-) of X is of the form 

oo 

0(0 = — [ I l{H(u,t) G ■}t r - 1 dtn n - 1 (du) (1) 

nn n J J 

s™- 1 o 

with 7 > and r G (0, oo). We refer to 7 as the intensity and to r as the distance exponent 
of the hyperplane process X. Clearly, is rotation invariant for all r > 0, but translation 
invariant only for r = 1. Therefore, since A is a Poisson process, X is always isotropic but 
stationary only for r = 1. 
The random polytope 

Z := P| H-, 

H&X 

is the zero cell of the hyperplane process A, where H~ denotes the (almost surely uniquely 
determined) closed half-space bounded by H which contains the origin. Clearly, Zq depends on 
7 and r; although this dependence is not made explicit by our notation. 
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For the distance exponent r = n the zero cell Zq is equal in distribution to the typical cell of 
a stationary Poisson-Voronoi tessellation (see [ID])- More detailed information on the topic of 
random tessellations is provided, e.g., in |34j . [30] and [36]. Poisson-Voronoi tessellations have 
been studied extensively in the literature; see, for instance, [7J, [13], [20], [21], [22], [23], [24J, 
[2"5] . [27] , [3T?] and [3E]. Stationary Poisson hyperplane tessellations have been considered, e.g., 
in PH, [U], [IS], [UJ, [IB], US], [35], |], HB, [32], H2J. Recently, also non-stationary Poisson 
hyperplane tessellations have attracked some attention (cf. [32], [ID])- A review of recent results 
on random polytopes is given in [3], [3]], [33], [3H Chapter 8]. 



3 A general formula for the variance 

Let X be a Poisson hyperplane process in W 1 with an intensity measure of the form 0. We 
assume that 7 > and r e (0, 00). By Fubini's theorem and basic properties of a Poisson 
process, we have 



E[V n (Z ) k ] =E I l Zo {x x ) dx x --- J l Zo (x k ) dx k 

(xt,...,x k e Z )dx 1 ...dx k 

00 

2 7 



exp 



S"- 1 
ir—lj-i un-li 



l{H(u,t) n [o,xx,...,x k ] ^ 0} 



x t^dtn^idu) 



dxi . . . dx k . 



(2) 



From Q we now deduce a lower and an upper estimate of the moments of V n (Zo). For e G S n 1 , 
we define 



c(n, r) 



r(r±2)' 



which is indeed independent of the choice of the unit vector e. The explicit value is determined 
by a suitable decomposition of spherical Lebesgue measure. 

The following result provides estimates from above and below for the moments of the volume 
of the zero cell. Note that the ratio of the upper and the lower bound is given by the ratio of 
the corresponding values of the Gamma functions in these estimates. 



Proposition 1. For k eN, we have 



n 



(. 



vi • 1 ) /■•:. . , 



>/.-! / P ( kfl , A .Jk ( _ UK n r 



In particular, for k = 1 we get 



Etv.(Zb)] = r(5 + i)K.(^ 



4 



Proof. Starting with (|2]), introducing polar coordinates and by symmetry, we get 



Since 



nv n (z ) k ] 

OO Si si 

= k 




exp 



27 



l{H(u, t) n [o, stvi, . . . , s k v k ] ^ 0} 



S"- 1 

x f^dtH^idu)] s*- 1 . . . s^ 1 U n ~\dv k ) . . . H n -\d Vl ) ds k ... d Sl 

00 Sl si 



< k 




exp 



27 



UKr. 



l{^(u,*)n[o,sn;i]^0} 



(S n - l ) k S"- 1 

tr— 1 jj. njn— 1 / J„,\ _ti— 1 „n— li/n— 1/ 



x f- 1 dtH n -\du) sr • ■ • s n k -"H n -\dv k ) . . . n n - L (dv!) ds k ... ds!. 



OO 

/ y l{if(u,*) n [o,siUi] 7^ %}t r - 1 dtn^idu) = ^c(n,r), 





the upper estimate follows easily. 

For the lower estimate, we observe that H(u, t) D [o, X\, ■ ■ ■ , x k ] 7^ if and only if there is 
some i G {1, • • • , k} such that H(u,t) D [o, Xj] 7^ 0. Hence, E[V^(Z ) fc ] can be estimated from 
below by 



exp 



2j_ 

7lK n 



J2H H (u,t) n [o,Xi] ^ 0K'- 1 dtU n -\du) 



dx\ . . . dx k 



exp 



S"- 1 

2tJNT 

nn„r 



c(n, r) 



Now the assertion can be shown by a straightforward calculation. 



□ 



Remark 1. The lower bound in Proposition[T]can also be obtained by an application of Holder's 
inequality. Moreover, for fixed r and k, the ratio of the upper and the lower bound is 

fkn \ , f n \ k . , fen 

r ( — + 1 j /T \- + lj = const(A;,r) • k~ . 

The following lemmas lead us successively to explicit formulas for the second moment and 
the variance of the volume of the zero cell. In a first step, the integral representation of 
the second moment of V n (Zo) from ^ will be simplified considerably by an application of a 
Blaschke-Petkantschin formula (cf. [31]). For this purpose, we define 



For 1 < m < n and Xi, . . . ,x m G M. n , let V m (xi, . . . , x m ) denote the m-dimensional volume 
of the parallelepiped spanned by these vectors. Further, denote by the unique rotation 
invariant probability measure on G(n,m). 
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Lemma 1. We have 

E[V n (Z ) 2 ] = h 



n:2 I exp 

span{ei,e2} 2 



2lL 

7lK n 



l{H(u,t) n [o,x 1 ,x 2 ] ^ 0}t 



r-l 



S 1 ™- 1 

xdm"- 1 (^)]v 2 (x 1 ,x 2 )"- 2 A s 2 pajl{eije2} (d(x 1 ,x 2 )). 

Proof. From (|2]), the linear Blaschke-Petkantschin formula [3U Theorem 7.2.1], the rotation 
invariance of spherical Lebesgue measure and the invariance of V 2 (-, •), we obtain that 



w n (z 



exp 



27 



UK,, 



l{H(u,t) D [o,x u x 2 ] ^ 0} 



G(n,2) L 2 S"- 1 

#r— 1 j + njn-1 



x t r - l dtU n -\du) V 2 (x u x 2 ) n ~ 2 X 2 L (d(x u x 2 )) v r 2 \dh) 

oo 

27 



''/-.2 / / eX P 

SO n span{ei,e2} 2 



1 1 h i 



i{H(u,t)n [0,^,^x2] ^0} 



5"" 1 

00 

27 



^.2 / exp 

SO n span{ei,e2} 2 



riK, n 



l{H(&U,t) D [0,&Fi,&F 2 ] 7^ 0} 





x f- 1 dtU n -\du)\ V 2 (x 1; x 2 ) n ~ 2 X 2 BpaD{eiM (d(x 1 , x 2 )) v n (dti) 

00 

27 



6„ i2 / exp 

span{ei,e 2 } 2 S" 1 " 1 

t.r—1 J-l njn— 1 / J„,\ V7 /„, „ \n-2 \2 



l{if(u,t)n[o,xi,x 2 ] ^0} 



which yields the assertion of the lemma. 



□ 



Next we simplify the inner double integral of the expression which was derived in Lemma [T] 
for E[V^(Z ) 2 ] by exploiting further the symmetry of the situation. 

Lemma 2. For x±,x 2 G span{ei,e 2 } C M n , we have 

oo 

J J l{H(u,t) n [o,x 1; x 2 ] ^ ^f^dtU^idu) 



S"- 1 



2-7T OO 



l< - 1 c(n,r) f [ l{H{(°°* 9 d ) ,t)n[o,x 1 ,x 2 ]^H}}t r - 1 dtd9. 








Proof. For n = 2 there is nothing to prove. Hence we can assume that n > 3. Let E :- 
span{ei,e 2 }. The map 



F : 



(mi,6 , ,u 2 ) i — y cos(6 i )mi + sin(6 l )u 2 , 



is injective and its image covers S n 1 up to a set of measure zero. Its Jacobian is 

JF(u u 6,u 2 ) = cos(fl)(sin(fl)) n - 3 , 
and hence the area-coarea formula (cf. p2, Theorem 3.2.22]) yields that 

oo 

J J l{H{u,t) n [o,x u x 2 ] ^ 0}f'" 1 dtH n -\du) 

7T 

2 OO 

= J J J J l{H{cos(e) Ul + sia{e)u 2 ,t) n [o,x 1 ,x 2 \ ^ 0} 
x cos{6){sm6) n - 3 t r - 1 dtU n - 3 {du 2 )d6U 1 {du 1 ). 

Since [o, ^1,^2] C E and 

t 



we get 



H(cos(6)u 1 + sin(0)« 2) t)fl£ = H E ui 



00 

y y i{H(u,t) n [o,x 1; x 2 ] ^ 0}r- i c/t^ i (rfM) 



COS6 1 



5™- 1 

7T 

2 OO 




(n-2)K n _ 2 / / / ) n[o,Xi,x 2 ] ^0} 



cos 6^ 

x cos(#) (sin 6) n -H r - l dt deU 1 {dm) 

7T 

2 OO 

n-2)K n _ 2 J y yi{fr B ( Ul ,*)n[o,zi,ar 2 ]^0} 
x (cos fl) r+1 (sin ^)"- 3 t r - 1 di d# ft 1 (dwi) 

7T 

2 

n-2)K n - 2 [ (cos 9) r+1 (sin 9) n - 3 d9 



sh 





OG 



y y l{#(wi, t) n [o,xi,x 2 ] ^ $}t r ~ 1 dtH 1 (du 1 ) 



x 

5 1 



2-7T OO 

^+J ) C (n,r) y y i{#( ,t) n [ 0>a;i) x 2 ] ^ ®}t^dtde, 





which completes the proof of the lemma. □ 

Having simplified the inner integral of the expression found in Lemma [I] for E[V^(Zq) 2 ], 
we reduce the outer integral in the next lemma by again taking advantage of the problem's 
symmetry. 



Lemma 3. We have 

E[V n (Z ) 2 ] =8nb n<2 



TT OO S 




exp 



2 7 c(n,r) r(| + l) 



riK r 



ooo 

2lT OO 



r+1 ' 
2 - 




l{ff( 



COS# 

sin 6 1 



, t) n [o, . 



COS <f\ I u 

sin^ J ' I 



] ^ 0}t r - 1 dtd0 



o o 

x s n ~ 1 u n ~ 1 (sintp) n ~ 2 dudsd<p. 
Proof. Combining Lemma [T] and Lemma [2j and introducing polar coordinates, we get 

nv n (z ) 2 } 

2ir oo oo 

2 7 c(n,r) r(| + l) 



2-7T 2-7T OO OO 



V2 




exp 





27T oo 




COS# 

sin# 



, t) n [o, , 



COS<y2 

sin 



it 



COS -0 

sin ip 







x us\ sin(<£> — ip)us\ n duds dip dip 

2n 2n oo s 



>n,2 




exp 



2 7 c(n,r) r(f + l) 



oo oo 

27T OO 



vrr( 



r+1' 
2 - 








cosy? 
sin ip 



u 



cos ?/> 
sin-?/' 



] ^tyXf-^dtdO 



x -us | sin(y? — ifj)us\ n duds dip dip 



where the symmetry in u and s is used to justify the second equality Hence we derive 



E[V n {Z f 



2-7T 2-7T OO S 



2b 



n.2 




exp 



oooo 

2ir oo 



2 7 c(n,r) + 



1 1 h > 




l\H( 



cos(# — -0) 



, t) n [o, . 







x ^-V- 1 ! 



2i m s 



4irb 



n,2 




sm(0 — ip) 
sin(y9 — ip) | n ~ 2 du ds dip dip 
2 7 c(n,r) r(§ + l) 



COs(<£ — "0)\ / U 

sin(^- V) i ' V0 



exp 



ooo 

2ir oo 




cos 6 



cos y \ 1 11 



simp J ' V 



o o 

.n—l „n— 1 



x w n x s n A |sin</?| n 2 duds dip 



TV OO S 



8tt6. 



n,2 




exp 







2 7 c(n,r) r(| + l) 



] ^ $Xf- l dtd6 




2-7T OO 

o o 

x s n ~ 1 u n '~ 1 (sin ip) n ~ 2 duds dip. 
For the last equation, we used that 



2-k CO s 




exp 



7T 



2 7 c(n,r) r(f + 1) 

nK n v^r(^) 



27T OO 




i^((l S f).')n[°. 



, sin # 



, it, — 1 „r7. — 1 1 ,„|n— 2 



COS (/9 \ /It 



sin <W ' V 



x w n i s n | sin 99 1 " duds dip 

-2-k 00 s 

2 7 c(n,r) r(| + l) 




exp 



-7T 

27T OO 

X 




/ cos 6 



cos(— ip)\ ( u 



x u n - x s n - x \s\YY{-p))\ n -'\-\) duds dip 
2 7 c(n,r) r(| + l) 



] ^(Af- 1 dtde 



7T OO S 




exp 



000 

27T OO 




i^((H).')n[ 



sin 6 1 

x u n ' L s n ' L (sin p) n ~ 2 duds dip 

2 7 c(n,r) r(| + l) 



cos ip \ u 
"- s l sinJ'U 



] ^0U r - 1 dtd0 





,n— 1 „n— 1 / 




exp 



000 

2tt 00 




sin # 



.n— 1 „n— 1 /_• ,.\n- 2 



COS <£> \ lU 



x u n s n (siny?)™ dudsdp 

2 1C (n,r) r(| + l) 



7T OO S 




exp 



000 

2ir 00 



iT(=±l) 






,n— 1 _n— 1/ 



COS <£> \ / U 

siny? J ' I 



«a.s 1 . ■ 1 . 1 1 ] ^ r^tj 



x -u™ 1 s n 1 (siny) n 2 dudsdip, 
which completes the proof. 



□ 



In the expression found in Lemma [3] the indicator function depends on 9,t,s,u and ip. Its 
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support can be determined explicitly. This is used in the proof of the following lemma where 
it is shown how the integration with respect to t in Lemma [3] can be carried out. 

Lemma 4. For u,sE (0, oo), u < s and ip G (0, tt), 

2tt oo 


= -u r [ (cos9Yd9 + -s r [ (cos 9Yd9, 
r J r J 

-f «( Vt *)- v 

where 

, . / Z — COS Lp\ / 7T 7T \ , r . 

2) := arctan I g I G (- - , - J , O, z) G (0, tt) x [0, lj. 

Proof. We have to determine the support of the indicator function 
[0,2vr] x [0,oo) -> {0,1}, 

™ - ^ (:::) (:::)• (;)^ 

First, note that 

is satisfied if and only if 

t G [0,w(cos#)+] U [0,s(cos(# - ip)) + ]. 
Hence we have to compare (ucos9) + and (scos(0 — <p)) + . We observe that 

/ - — cos p\ ( u \ 

u cos 9 = s cos(p — 9) 9 = arctan = a w, - , 

\ smcp J V s/ 

and therefore (|3]) is satisfied if and only if 

(M e ([-2' a (^' s)) X [Mcos0])U([a (p, - J ,p + - ] x [0, scos(0 - (p)]). 
Now the integral can be easily computed. □ 

In the following remark, we collect some facts which are helpful for a proper understanding 
of the formulas for the second moment and the variance of the volume of the zero cell which 
are finally presented in Theorem [I] 

Remark 2. (1) Let (p G (0, tt) and z G [0, 1]. Then 

~\ ~ V ~ 2 = a<y¥>, °' ) ~ ^ ~ a<y¥>, ^ = 2 ' 



since z h-> z) is increasing on [0, 1]. 
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(2) A simple calculation shows that 



[co&OYde 



r(l + i) 



(3) By the definition of k„ and 6 n;2 , we have 



n K n = Anb n ^ 2 J (sin ip) n dip. 
o 



In view of the next theorem and the consequences which will be derived from it, we provide 
a description of (E[V^(Z )]) as a double integral, i.e. 

(E[v n (z )]y 



IT 1 



, 2 87rfr n ,2 r / 2n\ ( nn n r 



r J \2jc(n,r) 




j.n—1 







(t r +l) 



- (sin <p)^dt dtp. (4) 



This follows from the special case k = 1 of Proposition [TJ Remark |2j (3), by the substitution 
z = (t r + l)" 1 , and by basic calculations involving Beta and Gamma functions. 

For the statement of the formulae for the second moment and the variance of V n (Zo), we 
introduce the auxiliary function 



IY- + 1 

K{t,<p) :- 2 



/ a{<p,t) 



\ 



f J {cos6) r d6+ J {cos6) r d6 

V -f <x(<p,t)-<p j 



for (t, if) E [0, 1] x (0, 7r). In the following, we will use that 

\ < F r (t, p)<t r + l< 2, (t, <p) E [0, 1] X (0, 7T), 
which is implied by Remark [2j (1) and (2). 



(5) 



Theorem 1. Let X be a Poisson hyperplane process in W 1 with an intensity measure of the 
form (jlj with intensity 7 > and distance exponent r E (0, 00). Then 



and 



E[V n (Z 



871b 



n.2 



r 



m (- 



7T 1 



nK n r 



r \ r J \2~fc(n,r) J J J FJt.tp) 

00 




t 



n-l 



^(sin<^) n 2 dtd<p 



Var[K(^o)]=^r^V W 



\ r J \27c(n, r) 

7T 1 



X 




_>;! 







F r (t,ip)- (t r + iy 



t"- 1 (sin <p) n - 2 dt dtp. 
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Proof. The formula for E[V^(Zo) 2 ] is implied by Lemma [3j Lemma [4] and straightforward cal- 
culations. The formula for Var[Ki(^o) 2 ] then follows from Q. □ 

Remark 3. In Appendix [A] numerical calculations of the variance using the numerical in- 
tegration functions of MATHEMATICA® are added. In small dimensions (n = 2, 3, 4) the 
numerical estimates of the variance for varying r are plotted for 7 = 1. In this case observe 
that by Proposition HI and Stirling's formula (see, e.g., j2]) we have lim E[V^(Z ) fc ] = /t^ for 

r >r\ 

— 21 

k e N and arbitrary fixed n > 2. In addition we study the choice 7 = 2 e(Tr) (^"(r "i 
which implies E[V^(Zo)] = 1 and lim E[Ki(Zo) fc ] = 1 for k G N and arbitrary fixed n > 2 by 

r— >oo 

Proposition [TJ Theorem |2j which will be proved in the following section, and Stirling's formula 
show that for both choices of 7 we obtain Var[V^(Z )] = O(^). 

On the other hand for specific choices of r (r = 1, r = 0.5n, r = n, r = 2n) the nu- 
merical estimates of the variance for varying dimension n are plotted for 7 = 1 and 7 = 
2 "(^) O-Xr The high- dimensional limiting behaviour of the moments and the vari- 

ance is studied in the following section. 



4 Variance estimates 

In the next theorem, estimates for Var[V^(Zo)] are provided. In these estimates two auxiliary 
quantities, D(n,r) and E(n,r), to be defined below, are involved. They are examined more 
closely in the following lemmas. For a better understanding of E(n,r), we introduce the 
auxiliary function M(v, r), v G [— 7r/2, 7r/2], for which an upper and a lower estimate is provided 
in Lemma [5] In Lemma [6] we establish estimates for E(n,r) which are based on Lemma [5] 

Definition 1. For v G [— f , f ] and r G (0, 00), we define 



and 

E(n, r) := 




2 



Since < tp/2 < tp — a(tp, t) < n/2, if tp G (0, n), we can rewrite E(n, r) in terms of the function 
M(-,r) as 



E(n, r) 



nT(f) 



vrr 



n-l- 



(sin 99) 



n-2 



+ r -1 M(> - t), r)] eft 
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Theorem 2. With these definitions, we have 

E(n,r)D(n,r) < Var[V n (Z Q )] < E(n,r) D(n,r) 4^ +1 . 
Proof. The mean value theorem and dol) yield that 



TT 1 

1 1 ,„\n-2 




F r (t,^)v (r + i) 

x 

7T 1 

n 1 



^ \ t n ~ ' (sin <p) n -'dt dtp 



> , 
r 2 



*// ( f + 1 - F ^ t '^) r ~ 1 ( sin ^ n ~ 2 ^^- 



o o 



Similarly, again by the mean value theorem and ^ we have 

^ Ir-^sin^)"- 2 ^^ 



TT 1 




F r (t,^ (r + i) 

x 

7T 1 

n 2n 

< 4-2^ 
r 




r + i ~p r {t, ^ e- 1 {sm ^ ) n ~ 2 dtd<f. 





Now the assertion follows easily. □ 

In the next lemma, estimates for the function M(-,r) from Definition [l] are provided. Note 
that for fixed v G (0,7r/2), the upper and the lower bound exhibit the same asymptotic be- 
haviour as r — > oo. Information about the asymptotic behaviour of M(-,r) is important for 
understanding the behaviour of the auxiliary quantity E(n,r). 

Lemma 5. For v G [0, |] and r G (0, oo), we have 

c±(cosvyr--v)<M(v,r)<C 1 fc ° S ^ +1 



y/r V2 / — ' — y/r sin(t>) 

for some constants c, C > which are independent of v and r. 

Proof. For the proof it is sufficient to consider the case v G (0, |). As lower bound for the 
cosine in (v , |), we use the line through cos(f) at v and through at |, which is the graph of 
the function 

J ' v- | V 2/ 

As upper bound for the cosine in (u, |), we use the tangent line at v, which is the graph of 

/ COS V \ 

9 ^ - sin(v) (9 -v J. 

V sin v / 

Now the assertion can be seen easily. □ 

The following lemma provides upper and lower bounds for the auxiliary quantity E(n, r) 
from Definition [I] appearing in the estimates for Var[V^(Zo)] in Theorem [2j 
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Lemma 6. (a) Let r G (0, oo) be fixed. Then 

c < E(n,r) < 2 
with a constant c > which depends on r but not on n. 
(b) For all r G (0, oo), we have 



with constants c, C > which are independent of r and n. 

Proof, (a) Let r G (0, oo) be fixed. Remark |2j (2), implies that < M(v,r) < 1 for all 
v G [— 7r/2, 7r/2]. Hence we get 

r(f) ?, 1 1\ n 

'n-l- 



E(n,r) < n ~ \ 2 , N / (sin^)™" 2 [ — + - 1 = — — + 1 < 2. 



On the other hand, we can estimate E(n,r) from below by 

2 1 

f) 



£?(n,r) > n I (sin^)"" 2 / e^M^t)^) dtdcp 



o o 



2 



T(f) f . M(f,r) 



-r(".,' ) J v " r n + r 



> (2(1 + r)Y l M r) =: c> 0. 



(b) We now consider the case of a parameter r G (0, oo) which is not necessarily fixed and may 
depend on n. This case requires a more careful analysis. Our plan is to split the integrals with 
respect to cp and t and to make use of Lemma |5} A suitable choice of the parameters t* G (0, 1) 
and if* G (0, |), which determine where we split the integrals, is provided later. However, we 
should point out that our choice will be independent of r, n. First, we have the decomposition 



E(n, r) 



n-2 



n— 1 ' 



+ 



nr(§; 



sin<£>) 
sin <f) 



o 

n-2 



[r +r - 1 M(o;(v3,t),r)+r- 1 M(v?-Q;(¥',t),r)] dt<fy> 
{t n+r - 1 M(a((p,t),r)+t n - 1 M(<f-a((p,t),r)} dt 



+ / {t n+r - 1 M{a{if,t),r)+t n ~ 1 M{if-a{<f,t),r)}dt dip. 



For the subsequent estimate, we use Lemma |5j < M(-, •) < 1, the fact that (tp,t) h-> a(y, i) 
is increasing in each variable separately and, in addition, that <p — a(<p,t) > f*/2 for ip > ip*, 
since a(<p,t) < <p/2. Hence we get 
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E(n, r) 

- Ar(^ 1 ) J 1 rJ \n + r nj r 
o 

+ f}}\, /(sinv^)- 2 ( —(t*) n+r + -(t*) n ) dtp 



2 

if* 



l 



< + 1 V (f - <p\ n - 2) + ((t*r +r + M(a(<p\nr)) 
+ (t*)" + M(f ,r) 

< C-j= max {cos(a(y9*, **)), cos ) } r + C max{sin(y9*), f } n , 

where we already used that our particular choice of £*,y?* will be independent of r,n. In fact, 
it seems to be difficult to choose t* and ip* in an optimal way as functions of r, n. However, in 
the special case r = n the following particular choice turns out to be optimal. In any case, we 
fix 

2 * • ( 2 \ 

t := —=, to := arcsm — — . 

For this choice we get 

sin(v?*) = -= « 0.8944, cos (£) « 0.8507, cos(a(p*, f )) = -= « 0.8944. 
v5 ~ v5 

Therefore our final result for the upper estimate of E(n,r) is 

E ( n,r)<c( 



Vr VV5/ VV5/ 

with some constant C > 0, independent of r, n. 

For the lower bound, we cut off the integral with respect to ip at ip° := 2 arcsin(;4=). Then 
we have 

sin(^°) = — « 0.9428, cos g)=^» 0.8165. 
3 V3 

It is not clear how to calculate the optimal choice of <p° for r depending on n in general. At 
least our choice is optimal for ip° e (0, |) and r = n. In any case, we thus obtain 

E(n, r) 

v>° i 

^ ^r ( L) J {s[nif)n ~ 2 J ^ +r " 1 ^(«(^t),0+^- 1 M(v^-a(^,t),r)] dtd V 
2 o o 

^° i 



> n J^L) / (sin ^ )ri ~ 2 / r+r_iM (f '0 ^ 
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n 



-M{l-y a ,n-2)M(^,r) 



n + r 

> c-^(sm<p°) n (cosf )' 



nr 



where we used that for < if < ip°, we have a(ip, t) < a(ip°, t) < and therefore M(a(ip, t),r) > 
M(^-, r). Thus we finally arrive at 

h(n,r) > c 



Rvr \ 3 / I v/3 

which concludes the proof of the lemma. □ 



In the following corollary, the estimates in ([6| are obtained from Proposition 
Stirling's formula. Similarly, the estimates in rt7l) are a consequence of Theorem 
and Stirling's formula. 



Hlb; 
q|2) 



by applying 
Lemma |6 



Corollary 1. For k G N, fixed r G (0, oo) and some constants c, C > 0, depending on r but 
not on n, we have 

kn kn 

c U(r)^ (l + £Tj <E[V n (Z ) k ]<Cn^ U(r)^ (l + (6) 

and 

A(r)^(l + ^) 5 V' < Var[V n (Z )} < C ^ U(r)^ (l + £) 5 V , (7) 



r + l 

7T 2 



eT(=±l)- 



A(r) 



For constant intensity 7, we infer from Corollary [l] that E[V^(Z ) fc ] and V^ar[V^(Z )] are 
divergent for fixed r and n — > 00. Hence, in the case of fixed distance exponent r G (0, 00) and 
constant intensity, the moments E[V n (Zo) fc ] are divergent. It is natural to ask if the intensity 7 
can be chosen as a function of n and r such that E[V„(Zo)] is equal to a positive constant A -1 . 
By Proposition [T] we have 

n \ I nn„r 



wy n {z Q )\ = r + 1) « 

Therefore, if we define 



270(71, r) 



with A > 0, then 

E[V n (Z )} = i. 

If we plug 7(r, n) into the lower estimate from Theorem [2j we deduce from Lemma [6] and 
Stirling's formula that (for fixed r) 

•p / 2n _|_ j\ ^ 

Var[K(Z )] > c ■ — ^ — — ^2 - ^ > c(r)A _2 v / n -> 00 as n -)■ 00. 



A 2 r(n + i)< 



16 



The preceding analysis suggests that in order to arrive at a limiting behaviour comparable 
to the case of a Poisson-Voronoi tessellation, i.e. with the variance converging to zero as n — > oo, 
we cannot choose the distance parameter r fixed but have to adjust it to the dimension. In 
fact, in the following theorem we consider the case where r is proportional to n which is the 
natural choice in view of the estimates that have been obtained. 

Corollary 2. Let r = an with constant a G (0, oo) and k G N. Then there are constants 
c,C > 0, independent ofn, such that 

kn kn 

Bm, <E[V,(Z,)']< C A»- ^" (8) 



k 



?y a \ Th J nj a \ 



and 




(9) 



B(a] 



n 



where 



Hence, if the intensity 7 is constant and r = an, thenE[V n (Z ) k ] andVax[V n (Z )] are converging 
to zero as n — >■ 00. 

Proof. The inequalities in ^ follow from Proposition [l] by means of Stirling's formula. 

The inequalities in ^ follow from Theorem |2j by estimating D(n, an) by means of Stirling's 
formula, and by estimating E(n, an) by means of Lemma M (b). □ 



Now we choose again r = an, with fixed a > 0, and determine the intensity 7 in such a way 
that the expected volume of the zero cell is equal to a positive constant. This is possible for a 
special intensity function depending on a, the dimension n and a positive constant A that can 
be prescribed arbitrarily. In the following theorem, we describe the asymptotic behaviour of 
the volume of the zero cell in such a setting. 



Theorem 3. Let r = an, with constant a G (0, 00) , and let the intensity be chosen as 

dTi^ K, f { 1 \ \ a 

70, n) = — — ^- Ar - + 1 ) n r 
2c(n,an) \ \a J 

with A > 0. Then the following is true. 

(a) For k G N, 

1 y(- + 1) 

<E[V n (Z ) k ] < 



\k - L nv uy j _ Afcr (i + 1 ) fcJ 

in particular, 

E[K(^o)] = \- 
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(b) There are constants c,C > 0, independent of n, such that 




< Var[V n (Z )] < C 



in particular, 



lim Var[V n (Z )] = 0. 



Proof. The inequalities in (a) follow from Lemma [T] and the special choice of the intensity 
as j(a,n). The inequalities in (b) follow from Theorem [2j since D(n,an) is constant for the 
intensity chosen as 7(0, n) and by estimating E(n, an) by means of Lemma M (b). □ 
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A Numerical Calculations 




50 100 150 200 250 

r 



Figure 1: Numerical evaluation of the formula for the variance from Theorem [T] using the 
numerical integration functions of MATHEMATICA®. In the top panel, we fix 7 = 1 and 
E[T^ l (Zo)] — lim E[Xi(Zo)] = E[V^(Z )] — n n is plotted as a solid line in the corresponding color 

r— »oo 

for comparison. In the bottom panel, 7 is chosen in such a way that E[V^(Z )] = 1. 
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Figure 2: Numerical evaluation of the formula for the variance from Theorem [T] using the 
numerical integration functions of MATHEMATICA®. In the top panel, we fix 7 = 1 and 
E[V^(Zq)] is plotted as a solid line in the corresponding color for comparison. In the bottom 
panel, 7 is chosen in such a way that E[V^(Z )] = 1. 
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